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ABSTRACT 



x 

Context. The unprecedented sensitivity of current X-ray telescopes allows for the first time to address the issue of the Fe K line 
complex variability patterns in bright, nearby AGNs. These kind of studies have the potential to map the accretion flow in the strong 
gravity regime of supermassive black holes. 

Aims. We examine XMM-Newton observations of the brightest sources of the FERO sample of radio-quiet type 1 AGNs (for a total 
^ . of 72 observations) with the aim of characterizing the temporal behaviour of Fe K complex features. 

Methods. A systematic mapping of residual flux above and below the continuum in the 4-9 keV range is performed in the time 
vs energy domain, with the purpose of identifying interesting spectral features in the three energy bands: 5.4-6.1 keV, 6.1-6.8 keV 
and 6.8-7.2 keV, corresponding respectively to the redshifted, rest frame and blueshifted or highly ionized Fe Ka line bands. The 
variability significance of rest frame and energy shifted Fe K lines is assessed by extracting light curves and comparing them with 
Monte Carlo simulations. 

^ ' Results. The time-averaged profile of the Fe K complex revealed spectral complexity in several observations. Red- and blue- shifted 

components (either in emission or absorption) were observed in 30 out of 72 observations, with an average (EW) ~ 90 eV for emission 
and (EW) ~ -30 eV for absorption features. We detected significant line variability (with confidence levels ranging between 90% and 
99.7%) within at least one of the above energy bands in 26 out of 72 observations on time scales of A^~6-30 ks. Reliability of these 
features has been carefully calculated using this sample and has been assessed at ~3cr confidence level. 

Conclusions. This work increases the currently scanty number of detections of variable, energy shifted, Fe lines and confirms the 
reliability of the claimed detections. We found that the distribution of detected features is peaked at high variability significances in 
the red- and blue- shifted energy bands rather than at rest-frame energies, suggesting an origin in a relativistically modified accretion 
flow. 
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^,1. Introduction NOGR07) demonstrated that Fe K line relativistic broadening 

is quite common, being effective in ~45% of the cases. 

Reflection features observed in the X-ray spectra of many Active In recent yearS; more diversity has been added t0 the general 

Galactic Nuclei (AGNs) represent a powerful tool to study the picture by the discovery of narrow and (apparently) transient 

properties of the matter flow accreting onto Supermassive Black f eatures j n the E-4-6 keV energy range (e.g. Turner et al. 2002; 

Holes (SMBH). In particular the fluorescent Fe K line can be Long i nott i et al. 2007a and references therein; Vaughan & Uttley 

used as a probe of the innermost regions of AGNs (e.g. Fabian et 2 008), characterized by relatively short-time scale variability, 

al 1989; Fabian et al. 2000; Reynolds & Nowak 2003; Miniutti of the order of tens of ks (e g Y aqoob et al. 2003; Iwasawa 

& Fabian 2005), where the effects of general and special rel- et al 2 004, hereafter IMF04; Tombesi et al. 2007, Petrucci et 

ativity become significant. Indeed, the relativistically skewed, al 2007) Xheir nature is not yet understood . Vaughan & Uttley 

redshifted and broadened profile of the line and its variability (2008) casted major doubts about the rel i ability of such shi f ted 

pattern strongly depend on the physical conditions of the accret- features , pointing out that many of the detections reported in 

ing matter and of the SMBH (e.g its spin). By now about 25% of the literature are not statistically significant and most probably 

all AGNs observed by XMM-Newton appear to show the pres- drawn from random fluctuations in the bulk of analysed X-ray 

ence of a relativistic Fe line in their X-ray spectra (Guainazzi data As sugge sted by these authors, the detailed and deep inves- 

et al. 2006; hereafter GBD06). Moreover, from a systematic tigation of a well defined sample of objectS5 as the one presented 

study of 37 very high S/N XMM-Newton observations of nearby herC; is necessary in order t0 achieve a firm conclusion on this 



(z < 0.05) radio quiet AGN, Nandra et al. 2007 (hereafter 



Send offprint requests to: demarco@sissa.it 



issue. 

The importance of such features relies on the fact that their en- 
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ergy is similar to that typical of relativistic Fe K line compo- 
nents, suggesting that we are observing the peak of a variable 
relativistic line profile produced in an extended region of the disk 
(in such scenario, the red tail would be too faint to be detectable 
above the continuum). Alternatively, intrinsically narrow emis- 
sion lines could be produced in discrete regions of inflows or 
outflows and redshifted as a consequence of relativistic mo- 
tions. Moreover, recent results point towards an origin in small 
emitting regions within the inner accretion disk (e.g. IMF04, 
Tombesi et al. 2007, Petrucci et al. 2007), as predicted by several 
models, as the "hot" orbiting spot model (Nayakshin & Kazanas 
2001; Dovciak et al. 2004) or the lamp-post model (Martocchia 
& Matt 1996; Dabrowski & Lasenby 2001; Miniutti & Fabian 
2004). Unfortunately, despite the unprecedented sensitivity of 
current X-ray telescopes, attempts to test these models are ham- 
pered by the too fine spectral features predicted (Goosmann et 
al. 2007). In this context, however, variability studies could be a 
key tool in providing insight onto their real nature (e.g. Ponti et 
al. 2004; IMF04; Miller et al. 2006; Tombesi et al. 2007). 
In this work we present the results obtained from the time- 
resolved X-ray spectral analysis of a well-defined sample (94% 
complete) of radio-quiet AGNs selected from the XMM-Newton 
archive. The sample is a subset of the FERO sample (FERO 
stands for "Finding Extreme Relativistic Objects") and was orig- 
inally selected by GBD06 to study the general properties of rela- 
tivistically broadened Fe Ka lines. We consider the same dataset 
with the aim of looking for transient features in the Fe K band 
and analyse their variability patterns by exploiting a uniform 
analysis of their highest quality observations. Several time re- 
solved studies have so far highlighted the variable nature of the 
Fe K line complex (e.g. Ponti et al. 2004; IMF04; De Marco et 
al. 2006; Miller et al. 2006; Petrucci et al. 2007; Tombesi et al. 
2007). Up to now, there is, however, no statistically solid charac- 
terization of its temporal behaviour. For the time-resolved anal- 
ysis we adopted the method based on the mapping of the excess 
residuals above and below the continuum in the time vs energy 
plane; this technique has been already widely used (Ballantyne 
et al. 2004; IMF04; Turner et al. 2006; Miller et al. 2006; Porquet 
et al. 2007; Tombesi et al. 2007; Petrucci et al. 2007) and proved 
to be very useful to characterize such variable X-ray spectral 
features. 

2. Selection of the sample 

The sample analysed in this paper was first presented in GBD06 
and comprises 33 poorly absorbed (namely with intrinsic N# < 
2 - 3 x 10 22 cm" 2 ) radio quiet AGNs selected from the RXTE 
Slew Survey (XSS, Revnivtsev et al. 2004) to have F2-\okev ^ 
1.5 x 10" 11 erg cm" 2 s _1 . To carry out a meaningful timing anal- 
ysis, all the XMM-Newton observations with EPIC pn exposure 
>10 ks were considered. 

This selection excluded the sources: HI 846-786 (being observed 
two times for a duration of <10 ks) and UGC 10683 (because no 
XMM-Newton observations were available as of 2009 February 
1). Furthermore, the source H0557-385 was discarded as its hard 
X-ray emission has been shown to be highly absorbed during the 
two >10 ks XMM-Newton observations (the spectral curvature 
below 6 keV implies the presence of an intervening neutral gas 
cloud, partially covering the X-ray source, with a column density 
Nh ~8xl0 23 cm~ 2 , Longinotti et al. 2009). Four observations 
(MCG-5-23-16, ID 0112830301; MRK 279, ID 0083960101; 
AKN 564, ID 0006810301; NGC 4593, ID 0109970101) were 
neglected after Good Time Interval (GTI) filtering, as they result 
badly affected by proton flares, reducing the EPIC pn effective 



exposure to <10 ks. Finally a 15 ks observations of AKN 564 (ID 
0006810101) was rejected, as the effective exposure, although 
>10 ks, is too short with respect to the time resolution needed 
for the timing analysis. 

Our final sample consists of 30 sources and 72 observations. 
Information about the sample is given in Table [2 
In the following each observation is denoted by the correspond- 
ing XMM-Newton revolution number. Whenever more than one 
observation is carried out during the same revolution, a letter in- 
dex indicates the chronological order. 

3. Analysis 

In order to characterize the variability pattern of the Fe K fea- 
tures using the so-called excess map technique (IMF04), it is 
important to identify a-priori the strongest lines present in the 
spectrum. The first step is thus to define the properties of the 
average spectrum, and then to analyse its temporal behaviour. 

3.1. Data reduction and time-averaged spectral analysis 

The analysis has been carried out using only the EPIC pn camera 
data because of its higher sensitivity (about a factor of two, e.g. 
Watson et al. 2001) in the Fe K band, with respect to the EPIC 
MOS. We used the XMM-SAS v. 8.0.0 software with CCF re- 
lease as of October 2006 for data reduction, and the Iheasoft v. 
6.0.3 package for data analysis. Time intervals free from high 
background events were selected by calculating the time series 
average count rate \i and standard deviation cr c during strong 
flares-free periods, in the E > 10 keV energy band and filter- 
ing out all the events with count rate exceeding the (arbitrarily 
fixed) threshold of \i + 3<x c . This choice allows to avoid large 
data gaps in the middle of the observation, which is a necessary 
condition for a timing analysis to be properly carried out. For 
each object, the source photons were extracted from a circular 
region of 45 arcsec radius, while the background ones were col- 
lected from adjacent source-free rectangular regions. Using the 
SAS task epatplot we verified the presence of significant pile-up; 
whenever some degree of pile-up was found it was minimized 
adopting an annular region for source counts extractiorQ. 
Average spectra were produced for each source and analysed us- 
ing the XSPEC v. 11.2.3 software package. The spectral fitting 
was performed in the 4-9 keV energy band, where Fe spectral 
signatures are expected, using simple models (i.e. a power law 
for the continuum and Gaussian components for emission and 
absorption lines). Some residual curvature of the spectrum pro- 
duced by warm absorber systems can yet be detected in this en- 
ergy band; in order to correct the shape of the continuum for it, 
a cold absorption component was added to the model (clearly 
the derived column densities, listed in Table [2 do not have any 
physical meaning). It is worth noticing that the presence of nar- 
row red- and blue- shifted absorption lines (Dadina et al. 2005; 
Reeves et al. 2004; Risaliti et al 2005; Cappi 2006 and refer- 
ences therein; Braito et al. 2007; Tombesi et al. 2009; Cappi et 
al. 2009) in the chosen energy band, did not influence signifi- 
cantly the determination of the continuum spectral shape. 
We modelled emission and absorption lines in the time-averaged 

1 Pile-up was revealed in the following observations: IC 4329a 
rev. 210 (carried out in "FullWindow" mode), MCG-5-23-16 rev. 
363 (carried out in "FullWindow" mode) and 1099 (carried out 
in "LargeWindow" mode), NGC 3227 rev. 1279 (carried out in 
"LargeWindow" mode) and ESO 141-G055 rev. 1435, 1436 and 1445 
(carried out in "FullWindow" mode). 
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spectra with simple Gaussian templates: a Gaussian was added 
to the model whenever we found evidence for it from the 99% 
confidence contours for intensity vs energy. The best fit models 
are reported in Table [3] 

3.2. Results from time-averaged spectra 

The "core" of the Fe Ka emission line is detected in all the ob- 
servations (see Table [3]). This line results usually relatively nar- 
row with Icr width of the order of the EPIC pn spectral resolu- 
tion at these energies (i.e. -0.13 keV) or less. A broader line at 
E-6.4 keV (<x -0.2-0.4 keV) is required in the spectral fitting 
of: NGC 3516 (rev. 1251), MCG -6-30-15 (rev. 108A, 301, 303), 
ARK 120 (rev. 679), MKN 766 (rev. 1000, 1001, 1003, 1004), 
MCG -2-58-22 (rev. 180), ESO 141-G055 (rev. 1435A, 1435B, 
1436), ESO 198-G024 (rev. 201), AKN 564 (rev. 930). The line 
energies are overall consistent with emission from neutral iron, 
although in some cases (NGC 3516 rev. 1251, 1253; MCG -6- 
30-15 (rev. 108 A); MKN 766 rev. 999, 1003; ESO 141-G055rev. 
1435B; ESO 198-G024 rev. 207; AKN 564 rev. 930) the implied 
energies might be indicative of mildly ionized matter (E>6.45 
keV) or more complex line profiles. 

Emission and/or absorption features at red- and blue- shifted en- 
ergies were observed in 30 out of 72 observations. We detected 
significant redshifted emission in 21 out of 72 observations (see 
Table with Icr width in the range <0.7 keV and (EW) ~ 90 
eV. It must be stressed that this is to be considered as a conser- 
vative number of detections because transient redshifted features 
might result smoothed out in the time averaged spectra analysis, 
and hence not significantly detectable. Time resolved surveys are 
the most suitable techniques for the study of this kind of features. 
Signatures of significant absorption lines blue ward the narrow 
"core" were instead detected in 15 out of 72 observations (see 
Table [B, at energies in the range E-6.6-8.5 keV, with (EW) ~ 
-30 eV. 

Finally, lines consistent with emission from highly ionized Fe 
Ka and/or Fe KJ3 were revealed in 28 out of 72 observations (see 
Table [T]). The Icr width of these lines were overall of the order 
of the energy resolution of the detector except for 12 observa- 
tions (MRK 509 rev. 250, 1168; IRAS 05078+1626 rev. 1410; 
MKN 766 rev. 82, 265, 999, 1001, 1002; NGC 526A rev. 647; 
NGC 4051 rev. 541; ESO 141-G055 rev. 1445, NGC 7213 rev. 
269), where a broader component was detected (<x~0.2-0.7 keV, 
EW-80-160 eV), which might be indicative of a blending, ei- 
ther between the ionized Fe Ka and neutral Fe KJ3 or with the 
blue peak of a relativistic component. 

3.3. Spectral variability 

The spectral variability has been investigated by systematically 
mapping residuals to the continuum emission in the time vs en- 
ergy domain over the 4-9 keV band. As most of the features ob- 
served in the time-averaged spectra are unresolved in their fine 
structure, the chosen energy resolution for the time-energy 'map- 
ping' is 0.1 keV for nearly all the source£| (i.e. approximately 
equal to the detector resolution in the analysed band). Though 
not producing a good oversampling for most of the features, this 
energy resolution represents the only viable choice for the avail- 
able data. The time resolution of the map has to be chosen in 
order to garantee a sufficient number of counts into each energy 



2 In the case of NGC 4051 (rev. 541), AKN 564 (rev. 930) and MRK 
704 (rev. 1074) we opted for a lower energy resolution (i.e. 0.2 keV, see 
Table [2]) in order to keep a high enough time resolution. 



bin. Hence, fixing the energy resolution, we extracted spectra at 
different time resolutions during the period of minimum flux of 
the source, and required a minimum of 20 counts per bin (for the 
X 2 statistics to be applicable) in the Fe K complex energy band 
(i.e. at least for E<7.5 keV). The adopted time resolutions are 
listed in Table [2] 

It is useful to compare the time resolutions used in the analysis 
with a typical time scale of the system. As indicative time-scale 
we can consider the estimated orbital period at 10 r g , given by 
(e.g. Bardeen et al. 1972): 

T orb = 3l0[a + (r/r g ) 3/2 ]M 7 (sec), 

r g being the gravitational radius, M 7 the mass of the black hole 
in units of 10 7 M o and a the adimensional angular momentum 
per unit mass. Following the Nyquist- Shannon sampling theo- 
rem, the condition At = T or b/2/3 (being ft >1 the sampling factor 
and At the chosen time resolution) must be satisfied for the cho- 
sen time resolution to oversample this characteristic time- scale. 
This is thus generally the case here that the estimated orbital 
time-scale at 10r^ is oversampled by the chosen time resolution. 
Only in four cases, being the orbital period very short (i.e. <2.5 
ks), the time resolution cannot adequately oversample it, as the 
minimum number of counts per energy bin must be preserved. 
The sources for which we were forced to undersample the orbital 
time-scale at 10r^ are: MCG-6-30-15, MRK 766, NGC 4051 and 
AKN 564. 

The method used to produce excess time-energy maps is 
extensively described in IMF04 and Tombesi et al. (2007) and 
briefly summarized in the following: 

i) each observation was divided into slices according to the cho- 

sen time resolution; 

ii) for each slice, the continuum was determined by excluding 
the Fe K energy band (i.e. typically at E-5-7 keV), and by 
rebinning the spectrum to have at least 50 counts per energy 
bin in order to garantee a reliable continuum determination 
in the high energy part of the spectrum (E>7.5 keV); 

iii) the baseline spectral model for the continuum is a power law 
plus cold absorption (with the column density N# parameter 
fixed at the value obtained from the average spectrum); 

iv) data were then rebinned at the energy resolution of the map. 
Residuals to the best fit continuum model were computed in 
the 4-9 keV energy band and visualized on a time vs energy 
plane; 

v) a circular Gaussian low-pass filter with cri p = 0.85 pixel is 
used to smooth the image in the time-energy plane so as to 
reduce random noise between adjacent pixels. 

In the maps, signal-to-noise (S/N) ratios of residual flux are dis- 
played (see Figs. H]-©. The noise includes contributions from 
uncertainties in the subtracted background and power law model 
components, and in the excess residual counts. 
The excess map technique is useful to identify spectrally 
and temporally transient, narrow features and trace their en- 
ergy/intensity temporal evolution. As these features could appear 
at an energy which is not known apriori we defined three physi- 
cally motivated energy bands for the extraction of residuals light 
curves, which allowed to infer the statistical significance of such 
structures in a uniform and systematic way: 

A: 5.4-6.1 keV, red shifted Fe Ka line band; 

B: 6.1-6.8 keV, neutral and/or mildly ionized Fe Ka line band; 

C: 6.8-7.2 keV, Fe KJ3 and/or highly ionized Fe Ka band. 
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Table 1. Time-averaged spectral properties of the observations. The revolution number of the observations where the most relevant 
features in the 4-9 keV spectra were detected is reported, including a flag indicating significances obtained from intensity vs energy 
contour plots. Multiple flags refer to the corresponding features in order of increasing energy. 



Source 


Redshifted Emission features 


Broad Fe Kd 


Highly Ionized Fe Kd/Fe KJ3 


Absorption features 


IC 4329a 


670* 




670* 


670** 


MCG -5-23-16 


363**, 1099* 




1099* 




NGC 3783 


193**, 371**, 372* 




193*, 371*, 372* 


193*'**, 372* 


NGC 3516 


245*, 352*, 1250** 


1251* 




1250*, 1251*, 1252*, 1253* 


A A TT> A/ C TXC\ 

MRK 509 






250 , 1073 , 1074 , 1168 


1073 , 1168 ' 


MCG -6-30-15 


108A**, 302**, 303* 


108A**, 301*, 303* 


303** 


108 A* 


a itr\r~^ i o 1 1 11 

MCG +8-11-11 






794** 




NGC 7314 






256* 




ARK 120 




679* 


679* 




MRK 279 








1088**'* 


NGC 3227 


1279* 








IRAS 05078+1626 






1410** 




MRK 590 






837* 




MRK 766 


265**, 999*, 1000**, 1001**, 1002**, 1003** 


1000*, 1001*, 1003*, 1004** 


82*, 265*, 999**, 1000*, 1001*, 1002* 




NGC 7469 






912** 




MCG -2-58-22 




180** 






NGC 526A 






647* 


647** 


NGC 4051 


263*, 541* 




541** 


263** 


ESO 141-G055 




1435A*, 1435B*, 1436* 


1436**, 1445* 


1445* 


ESO 198-G024 




207** 




1128** 


AKN 564 




930* 






NGC 7213 






269** 




NGC 4593 






465* 





(*): detection significance > Act; 

(**): detection significance between ~ 3cr-4<x; 



Despite the widths of the chosen energy bands are greater than 
the spectral binning of the maps, only relatively narrow features 
are picked up in the residuals analysis. Any other broad feature 
spreading across the entire width of each band would be ac- 
counted for by the continuum modelling for each time resolved 
spectrum. Moreover, adopting wider energy intervals allows to 
increase the statistics in the integrated residual counts. 
Light curves of the resulting excess flux and relative excess maps 
are reported for the observations where significant variability has 
been detected (see Figs.[4]-[9]). The light curves are renormalized 
for their corresponding average value. In the figures the maps are 
all shown at a time resolution of 2.5 ks, to better display excess 
residuals. The procedure used to calculate the light curves errors 
is outlined in Sect. 13.41 



with a Poisson noise distribution, which approaches the normal 
distribution whenever the number of counts per bin is >20. 
The square root of the simulated light curves mean variance (in 
the band of interest) is thus taken as error of the real light curve. 
This is justified by the fact that in the simulations every Gaussian 
component flux is assumed constant (and consequently equal to 
zero in case of no line detection), hence a measure of the vari- 
ance is a good estimate of the real time series errors. 
The variability significance has been calculated by comparing 
the real light curve variance (cr 2 .) with the simulated ones (cr 2 s ). If 
N is the number of simulated light curves for which cr 2 s >cr 2 , the 
variability significance is given by (1-N/N^ m ). Table U contains 
a list of the computed variability significances in the considered 
energy bands. 



3.4. Light curves error estimation and variability significance 

The excess map technique is useful to reveal the presence of nar- 
row transient features (IMF04) and is not affected by the pres- 
ence of broad line components (as the adopted continuum model 
would mimic the shape of such lines). The smoothing process 
applied (see previous section) prevents, however, the direct de- 
termination of the lines significance. For instance, the errors on 
the flux measurements in the time series cannot be estimated us- 
ing a counting statistics as, after the smoothing, adjacent pix- 
els in the map are no more independent. Hence we followed 
the procedure described in IMF04 to assess the errors via sim- 
ulations. We simulated N^ m =1000 time-energy maps, assum- 
ing constant spectral components (whose parameters are those 
of the average spectrum best-fit model, as reported in Table [3]) 
and apower law normalization varying as the 0.3-10 keV light 
curvfl These simulations are carried out using the fakeit Xspec 
task, which allows to create fake source and background spectra 



3 In the case of MCG-6-30-15 the curvature induced by the broad Fe 
line is treated as due to continuum emission (see Appendix). Since the 
broad line (E~3-7 keV) is less variable than the broad band 0.3-10 keV 



3.5. Robustness against continuum modelling and 
background filtering 

The technique just summarized assumes a simple modelization 
for the underlying continuum (i.e. a power law plus cold ab- 
sorption). As this assumption may affect the flux estimation of 
the variable features, we tested its robustness by considering 
a more complex continuum model, comprising an ionized ab- 
sorber component (absori model in Xspec, Done et al. 1992) 
rather than the cold absorption one. The test was performed for 
IC 4329a, the brightest source of the sample, with the highest 
statistics to provide the most accurate estimate of the contin- 
uum. For each energy band, we found that the calculated resid- 
uals are within the flux errors derived by the simple power law 
plus cold absorption continuum fit. On average, the difference 
between residuals calculated with the two continuum models is 
around 30% of the light curve error. This is due to several fac- 
tors: the narrow band considered for the analysis (E=4-9 keV), 
the high energy resolution of the maps (suitable for the study of 



continuum (e.g. Vaughan & Fabian 2004), the power law normalization 
is assumed to follow the 4-9 keV light curve in the simulated maps. 



B. De Marco et al.: Probing variability patterns of the Fe K line complex in bright nearby AGNs 



5 




-4 -2 2 4 6 
Fixed N (10" 5 ph/s/cm 2 ) 



















► ♦ 












Core : 



-2 2 4 6 
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Fig. 1. Comparison between line flux measurements for sub- 
tracted continuum model including a power law plus either a 
fixed or variable cold absorption N# parameter. The fluxes are 
obtained from residuals in the A, B and C energy bands (from 
the top to the bottom panel, respectively). 



narrow spectral features) and the lower variability level of sec- 
ondary components (i.e. reflection continuum and ionized ab- 
sorber systems) with respect to the primary power law. Hence, 
we conclude that the underlying continuum in the chosen band- 
pass is well represented by a simple power law, with cold ab- 
sorption being a good enough approximation for the modelling 
of the residual curvature. 

It is worth noting that the reflection continuum associated to the 
narrow and constant Fe Ka line dominates over that associated 
to any transient feature. This dominant reflection component will 
provide a constant contribution in the flux measurements. Hence 
it will not introduce spurious variability in the residuals maps. 
As a consequence of the adopted technique and high energy res- 
olution (i.e. AE =0.1-0.2 keV), the maps will be sensitive only 
to narrow spectral features (i.e. both spectral lines and narrow 
structures of the relativistic line profile.) We further checked if 
the assumption of a fixed value for the cold absorption N# pa- 
rameter in the continuum subtraction process (see Sect. 13.31) can 
produce spurious variable features (most of all in the redshifted 
energy band). This might be the case, for instance, in the pres- 
ence of any residual curvature induced by variable and/or com- 
plex absorbers (i.e. warm absorber) and/or more complex con- 



tinuum modelling. The test is exploited for the long-look (rev. 
301,302 and 303) of MCG -6-30-15, this source being charac- 
terized by a strongly variable X-ray flux and a complex spec- 
trum (including both a warm absorber and a relativistic Fe line). 
In FigsQ]we show a comparison between the residual line flux 
measurements (as obtained from the excess map analysis in band 
A, B and C) for the subtracted continuum model including (1) a 
variable or (2) a fixed N#. As expected, the line fluxes are well 
correlated (the linear Pearson correlation coefficients are 0.68, 
0.85 and 0.91 respectively for the three bands), although the scat- 
tering from the best fit linear model increases for lower energy 
bands, where the effects of continuum modelling are stronger. 
However the relative error due to the deviation in the line flux 
(l)-flux (2) plot from the linear model is smaller than the sta- 
tistical error in the line flux measurement as obtained from sim- 
ulations (see Sect. 13.4b . It is worth noticing that the scattering 
depends also on the energy and time resolutions adopted, de- 
creasing for lower resolutions. 

We conclude that our assumption of constant N# does not invali- 
date our results. Being the case of MCG -6-30-15 one of the most 
representative of a strongly variable bright source with broad rel- 
ativistic Fe line emission, we infer that this assumption can be 
considered valid for all the other sources of the sample. 

As some AGN show a spectral steepening when they get 
brighter, we also tested whether significant variations in the 
power law T parameter, which might influence the line flux mea- 
surements, are registered during each observation where signifi- 
cant variability in at least one of the three energy bands has been 
observed. The fractional deviation of the power law index from a 
constant model is always <20%. We rule out the possibility that 
such small deviations can produce spurious variations in the line 
fluxes, being the statistical errors on the fluxes generally larger. 
This result validates again our assumption of a constant T param- 
eter in the simulations (see Sect. 13.4b . Finally we exploited an 
a-posteriori check on the filtered background in order to rule out 
this can lead to fake results. Indeed energy-dependence of soft- 
p + flares might in principle produce spurious variable features in 
the spectrum of the source. To check this we inspected the back- 
ground 4-9 keV light curves during observations where signifi- 
cant variability in the Fe K band was found. After soft-p + flares 
removal, the background contributes less than 8% to the 4-9 keV 
count rate in most observations. Such a low background fraction 
is not expected to have any effect on our results. However dur- 
ing 7 of these observations the fraction of background count rate, 
even after soft-p + flares removal, remains high, falling between 
11% and 38% of the source one. To check if this may induce 
spurious energy dependent variability in the source spectrum we 
inspected the background spectra during periods of maximum 
count rate and compared them to the time averaged one. We fit- 
ted these spectra using a broken power law model, obtaining in 
all the cases a good fit. We found that the background spectral 
shape does not change, apart from a significant increase in nor- 
malization. This means that no energy dependence is present in 
flares. 



4. Results from excess maps 

The excess map technique revealed variability at >90% confi- 
dence level (against random noise) within at least one of the 
three (red, rest- frame and blue) energy bands in 26 out of 72 ob- 
servations (see Table |4]). Whereas, the overall number of energy 
bands showing excess variability in the sample is 34 with signif- 
icances in the range 90-99.7% each. In order to assess the sta- 
tistical significance of this result we carried out a binomial test. 
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Fig. 2. Distribution of variability significances in the three anal- 
ysed energy bands of the FOM- selected sub-sample (continu- 
ous curve). The dashed curve represents the distribution obtained 
from simulations whose variability is produced by random fluc- 
tuations. 
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Fig. 3. Distribution of variability significances of the FOM- 
selected sub- sample: comparison between the "core" (B band, 
dashed line) and the "red+blue" (bands A and C, continuous 
line) distributions. 



As we looked for variability within three distinguished energy 
bands (see Sect. 13.31) the total number of trials that must be taken 
into account in the computation is 216 (i.e. number of observa- 
tions times number of analysed energy bands). Furthermore, we 
make the conservative assumption that all the observed variable 
features are detected at the lower limit confidence threshold, i.e. 
90%. The test yielded 99.5%, i.e. about 3cr, as the probability 
for rejection of the null hypothesis. In other words, the num- 
ber of variable features detections in the sample is significantly 
greater than the one expected from a random distribution. In the 
former calculation we conservatively assumed the three energy 
bands to be independent. This is justified by the fact that, during 
the same observation, significant variability is rarely observed in 
more than one band (i.e. only in 7 out of 34 cases). Although 
from a physical point of view we would expect the three bands 
to be somewhat linked in their variability properties, we must 
indeed take into account several observational effects which al- 
low to consider the three bands independently. For example, the 
presence of blended constant components (this effect should be 
dominant in Band B, due to the presence of the "core" of the Fe 



Ka line, which plausibly arises in the outer regions of the disk 
or even farther), and the fact that S/N ratio decreases at higher 
energies, may result in an underestimate of the variability signif- 
icance in both B and C bands. 

As next step we tried to estimate and quantify potential observa- 
tional biases and/or limits that may hamper the detection of Fe K 
variability. One of the main observational issue in our analysis is 
the duration of the exposure. Physical considerations about the 
choice of the maps temporal sampling have been done referring 
to the orbital time scale at 10 r^ (See Sect. 13.31) . This is an arbi- 
trary choice, and does not address the possibility for variability 
phenomena to occur at shorter time scales. However, it repre- 
sents an observational bias which cannot be easily removed, as 
the quality of current data does not allow to sample shorter or- 
bital periods. In fact, in our analysis, sampling time scales must 
be chosen in order to preserve good statistics into each time re- 
solved spectrum, and should be able, at the same time, to probe 
short-term variability of the features in the Fe K band, even of 
relativistic ones. As relativistic effects are expected to be still 
observable at r~ 10 r g , the choice of considering the estimated 
orbital time scale at such distance as a characteristic time scale 
of the system is well justified. Unfortunately, for some of the 
analysed observations the total exposure is too short compared 
to the orbital time scale at 10 r^. This in principle forbids de- 
tection of features varying on time scales either of this order or 
longer. 

At the same time, the other important issue is the brightness of 
the source, as the possibility of obtaining significant detections 
depends on the total number of counts collected during the char- 
acteristic time scale of variations. 

In order to account for these two biases, we made a further se- 
lection on the sample. We defined the figure of merit (FOM): 

FOM = F2 - l0kevXT (erg cm- 2 ), 

T- orb 

as the total 2-10 keV flux per unit area during the whole ob- 
servation of duration T, normalized for the estimated orbital 
time scale T or t. This is equivalent to the total number of counts 
rescaled for the orbital time- scale. 

We arbitrarily sampled all the observations having 
FOM>l.4xlO~ 10 erg s" 1 cm -2 , which corresponds to the inclu- 
sion of those with an exposure 10 times longer than the orbital 
time scale and with a flux lower limit of F2_iofo?v>1.4xlO -11 
erg s" 1 cm -2 . This selection yields a total number of 33 "good" 
observation^ out of 72. Carrying out a binomial test on this 
sub-sample (which includes a total of 19 energy bands showing 
excess variability) we obtain a slightly higher probability, 
i.e. P=99.6%, that the recorded significant detections are 
not random. In the computation of FOM the main source of 
uncertainty is the BH mass, which affects the orbital time scale 
estimate. However, in most cases we used masses derived from 
optical/UV reverberation mapping measurements (see Ferrarese 
& Ford 2005 for a discussion on the reliability of reverberation 
mapping technique and a comparison with other methods). 
We collected different mass estimates (see references in Table 
0, obtained either via direct reverberation measurements or 
through secondary mass estimators based on reverberation map- 
ping, from recent literature (Wang & Zhang 2007, Peterson et al. 

4 IC 4329a (rev. 670), MCG -5-23-16 (rev. 1099), NGC 3783 (rev. 
193, 371, 372), NGC 3516 (rev. 245, 1250, 1251, 1252, 1253), MCG- 
6-30-15 (rev. 108A, 108B, 301, 302, 303), NGC 7314 (rev. 256, 1172), 
NGC 3227 (rev. 1279), MRK 766 (rev. 82, 265, 999, 1000, 1001, 1002, 
1003, 1004), NGC 7469 (rev. 912, 913), NGC 4051 (rev. 263, 541), 
MRK 110 (rev. 904), AKN 564 (rev. 930), NGC 4593 (rev. 465). 
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2004, Bian & Zhao 2003, Kaspi et al. 2000 and Ho 1998). For 
each source, we used the average mass value in the calculation 
of T orb (see Table©. In four cases (MCG -5-23-16, NGC 7314, 
MR 2251-178 and ESO 511-G030) reverberation mapping 
measurements were not available. Hence mass estimates from 
other methods were adopted (e.g. the photoionization method, 
Wandel & Mushotsky 1986, Padovani & Rafanelli 1988, or 
accretion disk models spectral fitting, Brunner et al. 2007). In 
those sources (MCG -6-30-15, NGC 526A, ESO 198-G024, 
AKN 564 and NGC 7213) where only one reverberation mass 
value was found in literature, available estimates from different 
methods were added in the computation of the average mass. 
Fig. [2] shows the distribution of variability significances in the 
three bands from this sub-sample. Assuming that statistical fluc- 
tuations dominate the observed variability in the data, we would 
expect the distribution of significances to be flat. Indeed, in this 
case, variances in the real data greater than those observed in the 
simulated light curves would be recorded 10% of the times with 
significances below 10%, 10% of the times with significances 
between 10%-20%, and so on, leading to a flat proability den- 
sity distribution. The dashed histogram in Fig. [2] illustrates the 
distribution of significances derived from simulated data whose 
variability is due to random fluctuations (assuming the same 
total number of observations as the real ones). We observe that 
the var-sig distribution is significantly peaked at high values, 
with a deficit at low significances, thus indicating the effective 
presence of the features we are catching. In particular, a K-S test 
strongly disfavours the observed distribution to be consistent 
with the random one (Pks =0.13%). Moreover, if we only 
focus on Band A and C where relativistic, energy shifted (either 
redwards or bluewards) features are expected to be detectable 
(Fig. [3j solid line), the clustering towards high significances 
results even more prominent (P^ rerf+fc/Me =0.06%). Conversely, the 
distribution is smoother and similar to the random one (Fig. O 
dashed line) when the energy band containing the narrow "core" 
of the E=6.4 keV Fe Ka component is isolated (P^ core =62%). 
This effect is indeed expected, since the narrow component 
contributes with an almost steady flux to the total band counts 
being most plausibly produced in regions far from the nucleus. 

5. Summary and discussion 

The main goal of this paper is to find and monitor transient emis- 
sion and absorption features in the Fe K line energy band of 
radio-quiet and bright AGNs. In this respect, as the detection of 
transient features is seriously biased by the lack of good statis- 
tics (time resolved spectra are usually characterized by low S/N, 
which prevents, in most cases to obtain high significance de- 
tections), the best method to derive strong conclusions on their 
reliability is by analysing a statistically complete sample. This 
issue is for the first time addressed in the present work, where 
we present both null and positive results from the uniform anal- 
ysis of an almost (94%) complete sample of sources, allowing to 
draw statistically meaningful conclusions. 
In this respect we studied a flux selected sub- sample of the 
FERO sample, defined by GBD06. The sample includes the 30 
brightest (in the RXTE Slew Survey, XSS, with a flux lower 
limit of F2-iofo?v > 1.5 x 10 _11 erg cm" 2 s _1 ) radio-quiet AGNs 
observed by XMM-Newton as of 1 February 2009. The observa- 
tions we analysed (for a total of 72) are chosen to have a duration 
>10 ks, in order to carry out a meaningful temporal analysis. 
Adopting the method of searching for excess variability (IMF04) 
in the energy bands of the redshifted (Band A), neutral-mildly 
ionized (Band B) and highly ionized Fe Ka and/or neutral Fe 



KJ3 (Band C) lines, we revealed significant (with a confidence 
level >90%) signatures of variable features in 26 out of 72 obser- 
vations, with probabilities in the range 90%-99.7%. Considering 
the total number of energy bands showing excess variability, this 
translates into 34 detections out 216 monitored energy bands. A 
binomial test for the significance of this result yielded a proba- 
bility of 99.5% (i.e. ~3<x) that the overall detections are not by 
chance. 

We further tested the detection frequency of variable features, 
after taking into account several observational biases (i.e. short 
duration of the exposure as compared with a typical dynamical 
time-scale of the system, as well as low flux observations) which 
might seriously prevent positive detections. We tried to remove 
these biases by defining a FOM ratio (see Sect. [4]). As these ob- 
servational limits forbid us to draw conclusions about Fe K vari- 
ability properties in the faintest sources of the sample, we will 
focus the following discussion on the FOM-selected dataset. 
Once the selection is accomplished, the distribution of variability 
significances in the overall three energy bands (with significant 
detections ranging from 91.7% to 99%), results peaked towards 
high values (i.e. >90%). Moreover, subtracting the contribution 
of the "core" energy band (band B), which is not expected to 
show signatures of strong variability due to the dominance of 
the constant Fe Ka emission, the skewness of the distribution 
towards high significances increases. All these results have been 
statistically checked through a K-S test, using a uniform distri- 
bution as the reference one. 

From the significances distribution analysis, it is reliable to as- 
sert that variability is commonly observed in the Fe K complex 
energy band of bright sources (see Fig. [2]), provided the typi- 
cal variability time scale is well sampled by the duration of the 
observation and the source is sufficiently bright. Infact, our esti- 
mated detection frequency is of 13 observations out of 33, show- 
ing variability in at least one of the three sampled energy bands, 
i.e. ~39%. Moreover, from Fig.[3]it is clear that the bulk of vari- 
ability comes from the red- and blue- shifted (with respect to 
the neutral "core") energy bands. This is most probably due to 
the fact that, if variability phenomena are still taking place in 
the neutral Fe Ka band, they result overwhelmed by the con- 
stancy (on short time scales) of the Fe line "core". Indeed, the 
detection frequency of variable features in the red- plus blue- 
shifted energy bands is 12 out of 33 observations, i.e. -36%. It 
is worth noting that this is very similar to the fraction of obser- 
vations (i.e. 12/37~32%) for which the characteristic emission 
radius of the relativistic component is constrained to be <50 r^ 
in the NOGR07 sample. 

Understanding the origin of the variable features in our sample 
from current data is not easy. Any excess of variability recorded 
in band A is most probably due to redshifted Fe K emission. 
Indeed this energy band, apart from the relativistic Fe K line 
components, is expected to include only the Fe Ka high order 
Compton shoulders, which, however, are predicted to be very 
faint and not observable with current instruments (e.g. George 
& Fabian 1991, Matt 2002). It has been claimed that partial cov- 
ering can induce spectral curvature in this band (e.g. Miller et 
al. 2008). Indeed, an ionized complex absorbing system might 
introduce further complexity and be responsible for variability 
as a consequence of variations in its properties (e.g. ionization 
state, covering factor, etc.). Although hard to predict, we expect 
that, in general, the partial covering gas would imprint its signa- 
ture in the B band as well as in the A and C bands. The plot of 
Fig. [3] conversely, clearly shows a lack of correlated variations 
between the "core" and the "red+blue" energy bands, disfavour- 
ing the partial covering interpretation. 
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On the other hand the other two bands, in addition to the con- 
stant (on short time scales) Fe Ka and Fe KJ3 emission compo- 
nents, can also include both the blue peak of a relativistic Fe line 
and ionized Fe K absorption features, which may be responsi- 
ble for the detected variability as well. Disentangling between 
the two cases would require a deeper analysis on any single 
case. However, residuals observed in band B are always posi- 
tive, which is proper of variable features in emission. On the 
contrary, there are three cases (NGC 3516 rev. 1250, 1251 and 
1252) where variable residuals in band C are negative, hence, 
most probably due to variable absorption features. This view is 
confirmed by the detection of two absorption lines in the time- 
averaged spectra of the three observations at these energies. In 
all the other cases the average residuals flux in band C is pos- 
itive, pointing to an origin in emitting material, most probably 
the peak of a relativistic Fe line. Assuming this is the viable in- 
terpretation, we can give 10/33, i.e. -30%, as an estimate for the 
detection frequency of relativistic variable features in bright and 
radio-quiet type 1 AGN, again of the same of order of the frac- 
tion deduced by NOGR07. 

Inspection of residuals light curves of the entire sample allows us 
to roughly compute the time scale at which the observed varia- 
tions take place. In almost all the cases these variations have not 
a regular shape. Only in three cases (IC 4329a rev. 670, NGC 
3783 rev. 372 and MRK 509 rev. 1073) they are characterized 
by a repeated, flares-like pattern. For these and other observa- 
tions, a deeper study can be addressed (i.e. NGC 3516, IMF04; 
MKN 766, Turner et al. 2006; NGC 3783, Tombesi et al. 2007; 
IC 4329a, De Marco et al. 2009), and we refer to the relative 
papers for an exhaustive discussion. In general, we decided to 
assume as the temporal time-scale of the variations, the time in- 
terval between the maximum peak in the light curve and either 
the following or the preceding flux minimum. The range of es- 
timated time- scales is A^=6-30 ks. This allowing to compute an 
upper limit on the extension of the region responsible for the 
observed variations. As variability from a source cannot be ob- 
served on time scales shorter than the light crossing time of the 
source itself, the upper limit on the size of the region is given 
by R < cAt (e.g. Mushotzky et al. 1993). If we consider the 
average (A^) -15.5 ks (which is also the average value of the 
FOM-selected sub-sample) as the reference one, the derived up- 
per limit on the region extension is R < 4.6 x 10 14 cm~140 for 
a black hole mass equal to the average value of the masses in the 
entire sample, i.e. {M BH ) ~ 2.24 x 10 7 M o . These calculations 
point to an origin into relatively small regions. 
Whether these regions are flares illuminated zones of the accre- 
tion disk or blobs of outflowing gas is not clear from our data, as 
overall the sensitivity of the observations does not allow to trace 
a clear pattern of variability. If associated to orbital motions, the 
observed variability would imply an average emitting radius of 
(r) ~7.9 r^ (assuming the average values of (At) -15.5 ks and 
(M B h) ~ 2.24 x 10 7 M o for the tim-scale and the black hole mass 
respectively). However, in most of the cases detection of signif- 
icant variability does not come with clear intensity modulations 
in the observed light curves, disfavouring the "hot spot" inter- 
pretation (e.g. Nayakshin & Kazanas 2001, Dovciak et al. 2004, 
Goosmann et al. 2007). Nevertheless, this does not completely 
rule out the model as the estimate of spot-induced modulations 
significance strongly depends on the possibility of observing the 
spot emission for a sufficiently long time (e.g. Vaughan 2005). 
This ultimately relies on the length of the observation and the 
permanence of the primary source flare. Moreover, the occur- 
rence of several flare events illuminating the disk is expected to 
generate more complex line patterns. 



While a discrete spot-like region of the disk would produce a 
smooth intensity/energy modulation pattern as a consequence of 
its orbital motion, this would not be the case if the emitting re- 
gion has a ring-like geometry. In that instance, the energy mod- 
ulations are expected to be weak, in agreement with results from 
extensive analysis of the variability patterns of some observa- 
tions of our sample (NGC 3783, Tombesi et al. 2007; IC 4329a, 
De Marco et al. 2009). Many authors pointed out the possibility 
to produce spiral density waves in non-self-gravitating models 
of accretion disks (e.g. Caunt & Tagger 2001; Hawley 2001). 
The effects of such structures on the Fe emission line would es- 
sentially be the production of several sub-peaks into the line pro- 
file, characterized by quasi-periodic variability as a consequence 
of the motion, and independent on continuum variations (e.g. 
Karas, Martocchia & Subr 2001, Hartnoll & Blackman 2002). 
Another suggested scenario is the production of line emission 
(via recombination of highly ionized Fe), during the ejection 
phase of blobs of gas originating from disk instabilities (Wang 
et al. 2000). In this case, detection of both red- and blue- shifted 
lines from both sides of the jet are expected. Some of the ob- 
served variable features in our sample are consistent with this 
picture, which is however not able to explain the modulations 
registered in some sources (e.g. NGC 3516, IMF04; MKN 766, 
Turner et al. 2006; NGC 3783, Tombesi et al. 2007; IC 4329a, 
De Marco et al. 2009). 

The present work is the first attempt to search for such patterns 
of variability. It makes use of the best quality data available up to 
now for this kind of studies. However, disentangling among pos- 
sible interpretations is still quite difficult. A qualitative step for- 
ward into the understanding of these issues will be provided by 
forthcoming missions, as the International X-ray Observatory 
(IXO). 
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Fig. 4. Left panels', data- to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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Fig. 5. Left panels', data-to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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Fig. 6. Left panels', data- to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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Fig. 7. Left panels', data-to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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Fig. 8. Left panels', data- to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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B. De Marco et al.: Probing variability patterns of the Fe K line complex in bright nearby AGNs 



Fig. 9. Left panels', data-to-model (power law plus cold absorption) ratios of the 4-9 keV time-averaged spectra; middle panels: S/N 
map of excess residuals in the time-energy plane (at a time resolution of 2500 s); right panels: 0.3-10 keV background subtracted 
continuum and residuals light curves (in bands A, B and C). The light curves are renormalized for the corresponding average flux. 
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Table 2. Sources list: (1) name; (2) RXTE count rate; (3) redshift; (4) estimated SMBH mass; (5) estimated Keplerian orbital period at r=10 v g \ 
(6) references for the mass estimates; (7) revolution number; (8) 2-10 keV flux; (9) duration of the exposure; (10) chosen time resolution for the 
excess map analysis; (11) chosen energy resolution for the excess map analysis; (12) FOM value (see Sect. |4j; (13) best fit value of the power 
law index; (14) best fit value of the cold absorption column density. The T and N H values are obtained assuming a simple "power law + cold 
absorption" for the fit of the time-averaged 4-9 keV continuum. The errors are computed at 90% confidence level. 



Sources 


RXTE 


z 


(Log(MWMo)) 


<T or ,(10r,)> 


Ref. 


Rev. Num. 


^2-10keV 


Obs. Duration 


At 


AE 


FOM 


r 


N„ 




(cts/s) 






(ks) 






(10 -11 erg cm -2 s" 1 ) 


(ks) 


(ks) 


(keV) 


(10 -11 erg cm -2 s" 1 ) 




(10 22 cm -2 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


IC 4329a 


7.29 


0.0161 


7.29 


19 


1,2,3,4 


210 


15.3 


10 


2.5 


0.1 


8.0 


2 - 04 -QQ8 


2 - 66+ 227 

35 +o:67 














670 


9.48 


133 


2.5 


0.1 


65.8 




MCG -5-23-16 


6.46 


0.0085 


7.60 


39 


13 


363 


7.43 


22 


2.5 


0.1 


4.2 


x 69 +o io 
i 70+8:0:3 


<3.50 














1099 


8.29 


122 


2.5 


0.1 


25.9 


2.67^2 


NGC 3783 


4.90 


0.0097 


7.01 


10 


1,2,3,4,5 


193 
371 


5.27 
3.96 


37 
131 


2.5 
2.5 


0.1 
0.1 


19.7 
52.3 


L85 -0.06 

1.67+0.20 


2.80^ 
1.61+1.10 














372 


5.56 


136 


2.5 


0.1 


76.3 


1 OO+0.03 
1 -° Z -0 07 


2 56 +0 - 50 

-0 70 


NGC 5548 


3.58 


0.0170 


7.71 


50 


1,2,3,4,5 


191 


3.14 


23 


2.5 


0.1 


1.4 


!- 67+ 006 

1.77+0.07 


<1.54 














290 


3.86 


93 


2.5 


0.1 


7.1 


<1.83 














291 


4.89 


37 


2.5 


0.1 


3.6 


1.89+0.10 


9 1 o+1.36 
A - LJ -\ 30 


NGC 3516 


3.21 


0.0088 


7.21 


16 


1,2,3,5 


245 


2.34 


129 


5.0 


0.1 


19.2 


1.81+0.09 


5.32 + \f 7 
7.82 + } :2 4 

4 - 76+ ?23 














352 


1.53 


128 


5.0 


0.1 


12.4 


1.82^-ji 














1250 


5.24 


52 


2.5 


0.1 


17.3 


2 15+0:06 














1251 


4.60 


69 


2.5 


0.1 


20.1 


9 9 1 +8:06 

z '- z ' i _Q 07 


6.07!^ 














1252 
1253 


3.64 
4.48 


68 
68 


2.5 
2.5 


0.1 
0.1 


15.7 
19.3 


9 90+0:09 

Z '- Z ' y -Q 08 

2 26 + ™ 7 


MRK 509 


3.12 


0.0344 


7.84 


68 


1,2,3,4,5 


161 


2.97 


30 


2.5 


0.1 


1.3 


1 63^ 

1 79+():i0 

L - ' A -0 15 

1.69+0.03 


<3.04 














250 


3.84 


43 


2.5 


0.1 


2.4 


<2.73 
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3.67 


85 


2.5 


0.1 


4.6 


<0.59 
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3.70 


46 


2.5 


0.1 


2.5 


1 74+0.09 

1 • ' ^_o 02 


<0.72 














1168 


4.26 


69 


2.5 


0.1 


4.4 


1 71+0:08 
A " -0.05 


<0.69 


MCG -6-30-15 


3.08 


0.0077 


5.84 


0.67 


3,8 


108A 


2.62 


43 


2.5 


0.1 


168.0 


2 01 + " W 
2.50+0.08 


<1.56 














108B 


3.60 


55 


2.5 


0.1 


295.3 


7.30+1.14 














301 


4.11 


84 


2.5 


0.1 


514.9 


2.47+0.05 


6.77+0.90 














302 


4.53 


127 


2.5 


0.1 


858.1 


IS 


±0.70 














303 


4.01 


125 


2.5 


0.1 


747.6 


9 1 7+u.zv 


L - J -0.87 


MCG +8-11-11 


2.36 


0.0205 


7.18 


15 


3 


794 


4.52 


38 


2.5 


0.1 


11.6 




<1.84 


NGC 7314 


2.16 


0.0048 


6.67 


4.6 


12 


256 


3.93 


43 


2.5 


0.1 


36.9 


2.12+0.09 


2.68 + |-g 














1172 


1.54 


82 


5.0 


0.1 


27.5 


1 94+0.14 
il-Al 


3 - 84 -14 


AKN 120 


2.14 


0.0323 


8.14 


134 


1,2,3,4,5 


679 


3.76 


112 


2.5 


0.1 


3.1 


1 07+O.O6 




MRK 279 


2.13 


0.0305 


7.79 


60 


1,2,3,5 


1087 


2.60 


59 


5.0 


0.1 


2.6 


1 91 +0.06 

-0.09 


<2.38 














1088 


2.43 


59 


5.0 


0.1 


2.4 


1 7 c +0.05 
A-/ -0.10 


<2.07 














1089 


2.38 


38 


5.0 


0.1 


1.5 


9 n^+o.08 

Z - UJ -0,15 


4 20 +2 06 

H - zu -2.00 


MR 2251-178 


2.10 


0.0640 


8.40 


246 


10 


447 


2.03 


64 


5.0 


0.1 


0.5 


1.61+0.09 


<2.14 


NGC 3227 


2.10 


0.0039 


7.39 


24 


1,2,3,4,5 


178 


0.85 


35 


5.0 


0.1 


1.2 


1 63 +U2U 

A ' DJ -0.13 


7 84 +2;/5 














1279 


3.64 


106 


2.5 


0.1 


16.1 


1.92+0.08 


2 89 +1 ° 8 


IRAS 05078+1626 


2.08 


0.0179 


7.04 


11 


1 


1410 


2.41 


58 


5.0 


0.1 


13.0 


1 7C+0.10 
i,/5 -0 09 


9 iro+1.46 
-1.25 


MRK 590 


1.95 


0.0264 


7.37 


23 


1,2,4,5 


837 


0.56 


108 


9.0 


0.1 


2.6 


2-30t^ 


c 9 1 +4.66 
JiZ,1 _l 97 


MRK 766 


1.77 


0.0129 


5.92 


0.82 


3 


82 


1.52 


39 


5.0 


0.1 


72.7 


9 J 

2 10 + fr°« 
L85|| 


3.46+2.29 














265 


2.44 


129 


5.0 


0.1 


386.0 


<1.20 














999 


0.72 


95 


8.0 


0.1 


83.9 


5 13+ 2 - 54 














1000 


1.13 


98 


6.5 


0.1 


135.8 


I 99+020 


<2.77 














1001 


1.40 


98 


6.0 


0.1 


168.3 


2 13+^ 
2 27+&1S 

2 01+^ 
2i +8:i 5 


<2.92 














1002 


1.74 


95 


5.0 


0.1 


202.7 


3 59 +L23 














1003 
1004 


1.48 
1.30 


98 
35 


5.5 
5.0 


0.1 
0.1 


177.9 
55.8 


<2.23 
<5.06 


NGC 7469 


1.69 


0.0164 


7.14 


13 


1,2,3,4,5 


192A 


2.49 


18 


5.0 


0.1 


3.4 


L88 -0.15 

1.79+0.11 


<3.97 














192B 


2.66 


23 


5.0 


0.1 


4.6 


<1.18 














912 


2.83 


85 


2.5 


0.1 


18.0 


1 99+8:8! 

L - yA -Q.Q7 


1 55 +118 
-0.93 














913 


2.93 


79 


2.5 


0.1 


17.3 


1 97+0.07 

— OS 


2.24+1.09 



Table 2. continued. 



Sources 
(1) 


RXTE 

(cts/s) 
(2) 


z 

(3) 


<Log(MWM )> 
(4) 


(ks) 
(5) 


Ref. 

(6) 


Rev. Num. 
(7) 


^2-lOkeV 

(10" 11 erg cm- 2 s" 1 ) 
(8) 


Obs. Duration 

(ks) 
(9) 


At 
(ks) 
(10) 


AE 
(keV) 
(11) 


FOM 
(10~ n erg cm -2 s _1 ) 
(12) 


r 

(13) 


N H 
(10 22 cm" 2 ) 
(14) 


MCG -2-58-22 


1.61 


0.0469 


8.00 


99 


1,3 


180 


3.02 


10 


2.5 


0.1 


0.3 


1.68:^ 


<4.04 


NGC 526A 


1.61 


0.0191 


6.98 


9 


1,13 


647 


2.20 


44 


2.5 


0.1 


10.5 


1 61 + ^ 


2-56!{i 


NGC 4051 


1.49 


0.0023 


5.90 


0.77 


2,3,4 


263 
541 


2.35 
0.56 


117 
50 


5.5 
6.0 


0.1 
0.2 


355.8 
36.2 


2 20+^' 


<4.35 


ESO 141-G055 


1.37 


0.0360 


7.51 


31 


1 


1435 A 
1435 B 
1436 
1445 


2.82 
2.54 
2.31 
2.88 


20 
21 
41 
78 


5.0 
5.0 
5.0 
5.0 


0.1 
0.1 
0.1 
0.1 


1.8 
1.7 

3.0 
7.2 


2 15 + ^ 
2.10 + ° : ^ 
i 94+0:06 

L y -0.08 


<3.88 
<5.42 
<3.31 
<0.78 


UGC 3973 


1.32 


0.0222 


7.74 


54 


1,2,4,5 


1247 
1263 
1332 
1535 


2.33 
2.01 
1.84 
0.77 


21 
20 
20 
87 


5.0 
5.0 
5.0 
6.0 


0.1 
0.1 
0.1 
0.1 


0.9 
0.7 
0.7 
1.2 


2 02 +Ui8 

-Q-17 

1 64+0- 16 
-Q.Q9 
1 63+ - 23 

1 70 + "- n 

-m 


<4.65 
<3.17 
<4.25 

5.87^32 


Fairall 9 


1.25 


0.0470 


8.00 


98 


1,2,3,4,5 


105 


1.16 


29 


5.0 


0.1 


0.3 


1 73 +uiy 


<5.00 


ESO 198-G024 


1.21 


0.0455 


8.38 


235 


1,9 


207 
1128 


1.30 
1.02 


31 
122 


5.0 
6.0 


0.1 
0.1 


0.2 
0.5 


1 84 +Uib 
i 7C+0.19 


<3.43 


MRK 110 


1.17 


0.0353 


6.99 


9.5 


2,3,4,5 


904 


2.80 


47 


2.5 


0.1 


13.8 


1.87±0.10 


1 7Q+1.38 


AKN 564 


1.13 


0.0247 


6.24 


1.7 


3,7 


930 


1.60 


99 


5.0 


0.2 


94.1 


2 50 +011 


<1.28 


NGC 7213 


1.12 


0.0058 


8.04 


107 


1,11 


269 


2.16 


42 


5.0 


0.1 


0.9 




<0.85 


ESO511-G030 


1.10 


0.0224 


8.66 


448 


6 


1402 


1.93 


109 


5.0 


0.1 


0.5 


1.89±0.08 


2 - 03 + )o7 


MRK 704 


1.06 


0.0292 


7.55 


35 




1074 
1630 


0.99 
1.07 


21 
98 


5.0 
5.5 


0.2 
0.1 


0.6 
3.0 


1.883$ 
1 07+0.19 


6.58^| 
<3.57 


NGC 4593 


1.05 


0.0047 


6.91 


7.9 


1,2,3,5 


465 


4.00 


76 


2.5 


0.1 


38.5 


1 84 +o:o6 


1 21 +0 - 91 



References: (1) Wang & Zhang 2007; (2) Peterson et al. 2004; (3) Bian & Zhao 2003; (4) Kaspi et al. 2000; (5) Ho 1998; (6) Winter et al. 2008; (7) Botte et al. 2004; (8) Vaughan et al. 2003; (9) Rokaki & Boisson 1999; (10) Brunner et al. 1997; 
(11) Nelson & Whittle 1995; (12) Padovani & Rafanelli 1988; (13) Wandel & Mushotzky 1986. 



Table 3. Gaussian components included in the best-fit models to the 4-9 keV time-averaged spectra, obtained assuming a "power law + cold 
absorption" as the continuum model. For each Gaussian component the line rest-frame centroid energy E (keV), width cr (keV), equivalent width 
EW (eV) and flux / (10~ 5 ph s" 1 cm -2 ) are reported. All the Gaussians are detected at >99% (A^ 2 >9.21) from confidence contours for line 
intensity vs energy. The errors are computed at 90% confidence level. 



Source 




Red 








Core 








Blue 








Absorption 






(rev. number) 


E 


cr 


EW 


/ 


E 


(T 


EW 


/ 


E 


(T 


EW 


/ 


E 


cr 


EW 


I 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(ID 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 




IC 4329a 
(rev. 210) 










6.43t£ 


<0.1 


CQ +24 
Jy -19 


9.9 


















IC 4329a 
(rev. 670) 


6 34 +010 


34 +0 - 17 


JZ -30 


5.9 


6.40±0.01 


0.06±0.03 


JZ -15 


5.7 


6 97+O.O6 


10 +008 


17 +17 
11 -9 


1.5 


7 67 +0 - 05 

/,D/ -0.02 


< 0.11 


"Hi 


-0.8 


MCG-5-23-16 

(rev. 363) 




°- 37 -a24 


74 -37 


7.6 


6 - 4z -am 


<0.1 




6.2 


















MCG-5-23-16 
(rev. 1099) 


6 35 +0 - 05 

"•■^ -0.06 


34 +0 - 07 

u -^-0.06 


90!}* 


9.6 


6.41±0.01 


<0.06 




4.4 


7.02+0.03 


<0.09 


13 -5 


1.1 










NGC 3783 
(rev. 193) 


6.21±0.05 


<0.13 


<28 


1.0 


6.38±0.01 


05 +U ii 


io9!j; 


6.8 


7.00+0.04 


<0.02 


26+9 


1.4 


6.60+0.04 
7.68+0.10 


<0.16 
0.20°_ 18 oq 


-16% 
-57+18 


-1.1 
-2.5 


NGC 3783 
(rev. 371) 


5 - 80 -a2o 


45 +u - 2? 

-0.18 


66- 


3.5 


6.39±0.02 


0.06±0.02 


137 -B 


6.4 


6.98+0.04 


0.10+0.04 


59+14 


2.3 










NGC 3783 
(rev. 372) 




46 +u - 65 

-0.05 




3.4 


6 40 +UU2 

«- i+u _0.02 


06 +u u4 


95! 9 4 


6.2 


7 01 +u -^ 


<0.16 


26- 


1.4 


6 67 +u - 25 

0,0 -0.19 


<0.47 


-20+7 


-1.4 


NGC 5548 










6.40±0.02 


<0.06 


/4 -16 


2.5 


















(rev. 191) 
































NGC5548 
(rev. 290) 










6.39±0.02 


O9 +0 03 

"•"^-0.02 


7 o+9 


3.3 


















NGC 5548 










6.37±0.03 


0.06< 0.11 


5V_\\ 


3.1 


















(rev. 291) 


































NGC 3516 
(rev. 245) 


6.29^ 


22 +aii 

u,zz -0.15 


03 -30 


2.3 


6.42±0.01 


<0.05 


88!- 


2.9 


















NGC 3516 
(rev. 352) 


6 40+o.ii 

"•^"-0.12 


62 +0 - 0K 

"•"-0.13 


214 +42 


6.3 


6.41±0.01 


0.07+0.01 


207+15 


5.2 


















NGC 3516 
(rev. 1250) 


5 98 +0 07 


< 0.17 


19+7 


1.4 


6.39±0.02 


O9 +0 03 

"•"^-0.02 


88±10 


5.5 










6.72+0.03 
7 3+o.04 

-003 


< 0.09 

< 0.09 


zz -8 

-27+8 


-1.3 
-1.3 


NGC 3516 
(rev. 1251) 










6 47 +0 07 

" ,H/ -0.03 


n i 0+0.07 
0.18_ 004 


221-f 


10.9 










6 66 +hl)1 

"•""-0.05 

7.04+0.06 


0.08+0.04 
< 0.08 


-69— 1 
-28!| 


-4.5 
-1.3 


NGC 3516 
(rev. 1252) 










6.39±0.01 


o.io-- 


134- 


6.8 










7.07*8? 


< 0.12 
<0.08 


-28+8 


-1.1 
-1.1 


NGC 3516 
(rev. 1253) 










6.46^ 


0- 14 -ao2 


184 +221 


9.1 










6.66 + |f 

7.08+S 


0.09+0.08 

0-015E 


_57+4V 7 
-27+10 


-3.4 
-1.2 


MRK 509 










6.38±0.03 


0.08±0.05 


'"-18 


2.6 


















(rev. 161) 


































MRK 509 
(rev. 250) 










6.42±0.03 


01 +a36 


30 +i/ 

-13 


1.4 


(. 64+0.26 
°- O4 -0.14 


23 +U i3 

w,zo -0.15 


48+24 


2.1 










MRK 509 
(rev. 1073) 










6.43±0.02 


0.12±0.03 


74±10eV 


2.8 


7 02 +U U6 

1 ,UZ -0.07 


<0.17 




0.8 


8.00-- 


<0.06 


-20— 
z "-n 


-0.5 


MRK 509 
(rev. 1074) 










-0.04 


12 +003 

" ,iZ -0.05 


76 +l2 


2.8 


6 9O +0 03 

"•^"-0.05 


<0.08 


27+10 


0.9 










MRK 509 
(rev. 1168) 










6.44±0.04 




39±12 


1.8 


6 91 +0 - 34 

°- yi -0.24 


66 +0 - 36 

"•""-0.39 


120!- 


4.7 


7 ^+0.04 

'• JJ -om 
8.50+0.05 


< 0.13 

< 0.12 


-26+J 1 

99 + 18 
P 


-1.0 
-0.6 


MCG-6-30-15 
(rev. 108A) 


5 84 +0 53 

J, ° -1.37 


U,J -0.47 




3.2 


6 45 +0 - 04 

D,HJ -0.10 


0.28-- 


162 +/3 


4.1 










6.67+0.04 


< 0.11 


_38 -205 


-1.1 


MCG-6-30-15 










6.38±0.02 


<0.09 


57±10 


2.5 


















(rev. 108B) 


































MCG-6-30-15 










6.43±0.04 


0.25±0.05 


iiZ -16 


5.3 


















(rev. 301) 
































MCG-6-30-15 

(rev. 302) 


5 84 +0 -^ 


21 +a28 
"• Zi -0.61 


25±10 


1.5 


6.42±0.03 


13 +00H 
"• -0.03 


67- 

"'-24 


3.3 


















MCG-6-30-15 

(rev. 303) 


c ^7+0.25 
J - D '-0.26 


071+ 0.22 
U - 7i -0.16 


146-f 


8.5 


6.40±0.03 


0.21±0.04 


130 +i6 


5.3 


6.90+0.04 


< 0.07 keV 


13- 


0.5 











Table 3. continued. 



Source 




Red 








Core 








Blue 








Absorption 






(rev. number) 


E 


cr 


EW 


/ 


E 


cr 


EW 


I 


E 


cr 


EW 


/ 


E 


cr 


EW 


/ 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) 


(17) 


MCG +8-11-11 

(rev. 794) 










6.42+0.01 


06 +U 02 

u ' uo -0.05 


119- 


5.9 


6.90321 


15 +(m 

u,iJ -0.73 


43+I8 


1.9 










NGC7314 
(rev. 256) 










6.40+0.03 


10 +UU5 


72+14 


3.1 


6.8633; 


12 +UU5 

U,iZ -0.06 


52+16 


1.9 










NGC7314 
(rev. 1172) 










6.42+0.02 


10 +OM 

u - -0.03 


134 +18 


2.4 


















ARK 120 
(rev. 679) 










6 41+ - 02 


15 +004 

u - 1J -0.07 


127 +y 

iz/ -14 


4.8 


6 94+ 014 

"•^-0.31 


12 +a26 

U - -0.05 


42+9 

HZ -10 


1.4 










MRK 279 










6.43+0.02 


0.08+0.03 


96 +i5 

yD -16 


2.6 


















(rev. 1087) 
































MRK 279 
(rev. 1088) 










6 41+0.O6 

°- Hi -0.07 


0.08tS 


09+23 


2.5 










6 69 +U 04 

-8-85 

7 76 +a07 


<0.08 
12 +015 

U - iZ -0.07 


-22 + * 
.55+18 


-0.6 
-1.0 


MRK 279 










6.41+0.03 


0.08+0.04 


113+21 


3.1 


















(rev. 1089) 


































MR 2251-178 
(rev. 447) 










6 42+ - 05 


0.10+0.08 


51 +14 


1.2 


















NGC 3227 










6.39+0.01 


0.05+0.02 




2.7 


















(rev. 178) 
































NGC 3227 
(rev. 1279) 


6 02 +U 08 

D,UZ -0.05 


<0.18 


27±9 


1.3 


6.40+0.01 


0.06+0.02 


iUZ -10 


4.2 


















IRAS 05078+1626 
(rev. 1410) 










^q+0.02 

°-^-0.03 


<0.10 


94+13 
* -33 


2.7 


6 66 +0 U9 

°- DD -0.32 


<0.35 


17+10 


0.5 










MKN 590 
(rev. 837) 










6.44+0.02 


U,U 1 -0.04 


i/z_ 31 


1.1 


7.02^ 


07 +a2 

U,U/ -0.02 


96+31 

yD -37 


0.5 










MKN 766 
(rev. 82) 










6 42 +u - U4 

U,H -0.03 


< 0.13 


50!g 


0.8 


6.80-- 


22 +u,1/ 

U,ZZ -0.12 


114^ 


1.6 










MKN 766 
(rev. 265) 


5 62 +0 - 04 

J,DZ -0.03 


< 0.71 




0.4 


f. oo+0.05 


< 0.10 


15!f 


0.4 


6 62 +0 31 

D,DZ -0.05 


32 +U U6 

U - JZ -0.14 


159- 


3.4 










MKN 766 
(rev. 999) 


5.98^ 


42 +u n 

U -^ Z -0.20 


185- 


1.9 


6.46+0.03 


< 0.13 


DZ -55 


0.7 


6 70 +u lu 

D,/U -0.13 


< 0.40 




0.4 










MKN 766 
(rev. 1000) 




°- 58 -ao9 


165^ 


2.2 


6.40+0.03 


<0.28 


H/ -24 


0.6 


6 62 +U U6 

°-° -0.09 


09 +U iU 

u - u ^-0.05 


56^ 


0.7 










MKN 766 
(rev. 1001) 




< 0.12 


19+10 


0.3 


6 41+0.15 

D,Hi -0.08 


< 0.28 


64+i« 

° -38 


0.9 


6 - 83 -ai9 


17 +0iy 
u - 17 -0.11 


"-25 


1.2 










MKN 766 
(rev. 1002) 


c 70+O.O6 


< 0.12 


18+9 


0.4 


6 40 +0 ° 2 


< 0.09 


OQ+14 

3 *-13 


0.7 


6 77+ U U6 

D,/ 7 -0.07 


20 +0 07 

U,ZU -0.06 


l° 6 -22 


1.7 










MKN 766 
(rev. 1003) 


c qo+0.03 
~ } - yj -0.02 


< 0.20 


20!- 


0.3 


6 54 +0 - 07 


0.26+0.06 


148+25 


2.0 


















MKN 766 
(rev. 1004) 










f. 44+0.20 

"•^-0.12 


27 +u - 12 

u,z/ -0.20 


139+115 


1.9 


















NGC 7469 










6.39+0.02 


<0.07 


108 -22 


2.8 


















(rev. 192A) 
































NGC 7469 










6.41+0.03 


< 0.14 


90- 


2.4 


















(rev. 192B) 
































NGC 7469 










6.40+0.01 


0.06+0.02 


95^ 


2.8 


7.00+0.05 


< 0.14 


21+9 


0.5 










(rev. 912) 


































NGC 7469 










6.42+0.02 


0.08+0.02 


94+12 
^-11 


2.9 


















(rev. 913) 
































MCG-2-58-22 
(rev. 180) 










6-39-S 


17+ Ui6 
u - 1 ' -0.08 


107^ 


3.6 


















NGC 526A 
(rev. 647) 










6.38^ 


<0.09 




1.0 


6.59tS 


32 +a22 

U - JZ -0.13 


82^ 


2.2 


7 97+0.04 
/,zz -0.05 


<0.12 


_ 19 +io 


-0.5 


NGC 4051 

(rev. 263) 


6 27 +0 - 24 

D,Z/ -0.20 


4 1+ 0.43 


68- 


1.9 


6.39+0.01 


<0.07 


68!- 


1.8 












<0.18 


ZH -12 


-0.5 


NGC 4051 
(rev. 541) 


6 13 +(m 


65 +016 

U,OJ -0.28 


485^ 


3.6 


6.44+0.02 


07 +O OJ 

u,u ' -0.04 


196 -36 


1.5 


7.08^ 


< 0.21 


49 -27 


0.3 











Table 3. continued. 



Source 




Red 






Core 








Blue 








Absorption 




(rev. number) 


E 


a EW 


/ 


E 


cr 


EW 


I 


E 


cr 


EW 


/ 


E 


cr 


EW / 


(1) 


U) 


(3) (4) 


(5) 


ip) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


(15) 


(16) (17) 


ESO 141-G055 
(rev. 1435 A) 








6 35 +(Ui 

U,J -0.09 


<0.27 


70 +3i 


2.0 
















ESO 141-G055 
(rev. 1435 B) 








6 54 +0 - 19 

U,J -0.20 


39 +0 - 2 ° 




5.1 
















ESO 141-G055 
(rev. 1436) 








6.35±0.07 


17 +008 

U - 11 -0.07 


1 17 +36 

11 '-33 


2.7 


6.90±0.03 


< 0.11 


5l + ' u 

-22 


1.0 








ESO 141-G055 
(rev. 1445) 








6 42 +0 - 03 

°-^ z -0.02 


< 0.09 


30+17 


0.9 


6 62 +0 - 2 ° 

°- DZ -0.17 


45 +u - 26 

U,HJ -0.33 


126 +35 

iZD -59 


3.5 


o ^i+0.03 


< 0.11 


-75±20 -1.3 



UGC 3973 
(rev. 1247) 



6.44±0.03 



<0.08 



69±20 



1.7 



°^-23 



pd 

b 

CD 

O 
O 
CD 



p 

I 

O 

Er 

CD 
CD 

CD 
O 

o 
3 



UGC 3973 
(rev. 1263) 



6.42±0.03 



AZZ -30 



UGC 3973 
(rev. 1332) 



6.38±0.04 0.11 ±0.04 



2.5 



284^ 



UGC 3973 
(rev. 1535) 



6.35±0.01 



0.11! 



3.1 



i ?6 +24 



Fairall 9 
(rev. 105) 



6.40±0.03 



< 0.11 



1.7 



ESO 198-G24 
(rev. 207) 



6.51+0.11 



<0.36 



60±29 



0.8 



07 +U U2 

U * -0.05 



_49+^ 

* y -25 



ESO 198-G024 
(rev 1128) 



6.41±0.02 



77 + 



0.9 



7.56±0.13 



<0.33 



-0.4 



MKN 110 
(rev. 904) 



6.41±0.04 



< 0.13 



46±13 



1.4 



y *-25 



AKN 564 
(rev. 930) 



6.54±0.10 



1.2 



6 74+0.13 



NGC 7213 
(rev. 269) 



6.40±0.02 



<0.08 



79+ 



59±24 



12 +004 

U,A -0.03 



ESO511-G030 
(rev. 1402) 



6.39±0.02 



95+13 



2.0 



6 37 +0 06 



13 +007 

U - i - , -0.06 



MRK 704 
(rev. 1074) 



159±43 



2.3 



MRK 704 
(rev. 1630) 



6.38±0.02 0.11 ±0.04 



140+ 



1.6 



6 93+0-05 



09 +U i2 



NGC 4593 
(rev 465) 



6.36±0.01 



< 0.11 



121+10 



5.2 



37+u 



CD 
2 



20 
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Appendix A: Notes on single sources 

The aim of this section is to discuss and compare our results 
with previous studies, as many of the sources of our sample have 
been extensively investigated in literature. In general, the best fit 
models we used as templates for the simulations (see Table [3} 
are in good agreement with the ones derived by other authors 
from deeper broad X-ray band analysis. We mostly refer to the 
NOGR07 survey for broad Fe lines, as they analysed many of 
the observations treated in this paper. In the following we will 
focus on the most relevant sources. 

IC 4329a\ The addition of a broad and redshifted emis- 
sion component to the best fit model is required at 99% confi- 
dence level (from intensity vs energy contour plots) during rev. 
670. The detection of this line is confirmed by previous stud- 
ies (Markowitz et al. 2006; NOGR07) on the same data set. 
However, its origin is controversial, being, most probably pro- 
duced at r>50 r g , where relativistic effects cannot be robustly 
measured (see NOGR07). We registered significant flux varia- 
tions within the energy band comprising this component, on time 
scales of ~32 ks. We refer to De Marco et al. 2009 for a more 
comprehensive study of such emission feature. 
On the contrary the redshifted component appears to be absent 
(in agreement with NOGR07 results) during the preceding ob- 
servation (rev. 210). In this case, the non-detection of the line 
may be a consequence of the short exposure (-10 ks). However, 
significant variability is detected in band B, pointing to the plau- 
sible presence of some kind of persisting, short- time scale vari- 
ability phenomenon taking place in this source. 
Detection of a blueshifted (~ 0.1c) and constant flux absorption 
line is registered during rev. 670, in agreement with Markowitz 
et al. 2006 and NOGR07, possibly ascribable to FeXXVI. 

MCG -5-23-16: This source has been deeply monitored via 
simultaneous observations with XMM-Newton, Chandra and 
Suzaku. This allowed accurate determination of the underlying 
continuum and detection of a relativistic Fe Ka line during rev. 
1099. We confirm detection of a broad component during both 
rev. 363 and 1099, after restricting the analysis energy band to 
E=4-9 keV. This line does not show significant variability, in 
good agreement with previous studies (Dewangan et al. 2003, 
Braito et al. 2007, NOGR07). 

NGC 3783: We revealed a redshifted emission line in all the 
three observations (rev. 193, 371, 372), in agreement with pub- 
lished results (Blustin et al. 2002; Reeves et al. 2004; NOGR07). 
However, NOGR07 associate the broad component observed 
during rev. 193 to emission originated at r>50r^. We also found 
signatures of a narrow absorption line at E-6.6 keV during 
rev. 103 and 372, ascribed by Reeves et al. 2004 to a blend of 
FeXXIII, FeXXIV and FeXXV lines at the systemic velocity of 
the source. Reeves et al. observed that the line decreased its in- 
tensity between rev. 371 and rev. 372, indicating a variability 
time scale of ~10 5 s for the highly ionized absorber. Although 
our technique does not test variability between different revolu- 
tions, we can nonetheless confirm these findings, as we do not 
have strong evidence for an absorption line detected during rev. 
371. However, it is worth noting that the line is present during 
rev. 193, in good agreement with results of NOGR07. 
An exhaustive temporal study of the spectral features in the 3 
XMM-Newton observations of NGC 3783 has been published 
by Tombesi et al. 2007. Our analysis resambled the overall re- 
sults of this paper. 

NGC 3516: Narrow, energy shifted emission components in 
the Fe K line profile were for the first time detected during a 
simultaneous Chandra and XMM-Newton (rev. 352) observa- 



tion of NGC 3516 (Turner et al. 2002). Variability of the nar- 
row redshifted emission feature at E-6.08 keV (Bianchi et al. 
2004) has been widely investigated by IMF04 for rev. 245 us- 
ing the excess map technique, finding a variability confidence 
level of 97%. The pre-defined redshifted energy band we used 
for our analysis (band A, see Sect. 13.3b is broader than the en- 
ergy band monitored by IMF04 (i.e. 5.8-6.2 keV), hence the light 
curve variations are smoothed out. However, restricting the en- 
ergy band to the one considered in IMF04 we achieve the same 
variability confidence level (i.e. 97.1%). 

In addition, different layers of absorbing gas were revealed dur- 
ing rev. 245 and 352 (Turner et al. 2005, NOGR07). Here we did 
not not account for these components as the narrow energy band 
in which we focus our analysis is not strongly affected by them. 
However, the technique we used to trace variability is nonethe- 
less able to reveal high energy transient components of the ab- 
sorbing gas. On the contrary, we observed two deep absorption 
lines in the average spectra of NGC 3516 during four consecu- 
tive revolutions (i.e. rev. 1250, 1251, 1252, 1253, covering the 
period 2006 Oct. 6-13), in the 6.7-7.1 keV regime (coincident 
with the C band), which can be associated to H-like and He-like 
species of Fe. These lines were ascribed by Turner et al. 2008 to 
a disk wind with outflowing velocities in the range 1000-2000 
km s" 1 . All these absorbing structures seem responsible for the 
significant variability revealed in band C during rev. 245, 1250, 
1251 and 1252. We point out that the lack of variability during 
rev. 1253 is not due to statistical problems as the observation has 
approximately the same duration and flux as the previous three. 

MRK509: As pointed out by NOGR07, the statistical signif- 
icance of the relativistic Fe Ka line during rev. 161 is compa- 
rable to the one obtained using a distant reflector model. Hence 
we adopted a simple, relatively narrow Gaussian template to fit 
the line. We further observed two absorption features during rev. 
1073 and 1168. One of these (E-7.33 keV) is also revealed in 
the synthesis spectrum (averaged out over all the five exposures) 
of Ponti et al. 2009, and its origin is ascribed to highly ion- 
ized outflowing (v~ 14000 km s" 1 ) gas. On the contrary the high 
energy one (E-8.0 keV) is also detected by Cappi et al. 2009 
and Tombesi et al. 2009. We observed significant variability in 
band B during rev. 161 and 1168, in agreement with results ob- 
tained from the total RMS spectrum analysis (Ponti et al. 2009). 
Moreover we found signatures of significant variability in Band 
A during rev. 161 and 1073. 

MCG -6-30-15: In this particular case the relativistic Fe Ka 
line (e.g. Wilms et al. 2001) extends over the entire energy band 
where we focused our analysis. As the band is too narrow to ob- 
tain a good fit of the broad Fe line, the induced residual spec- 
tral bending can be modelled as due to continuum emission. 
Hence, as a consequence of the adopted technique and energy 
resolution (i.e. AE =0.1 keV), the mapping procedure is sensi- 
tive only to narrow spectral features above the broad line, and 
the corresponding flux measurements are robust against contin- 
uum modelling. The RMS spectrum (with a time resolution of 6 
ks) of the two observations carried out during rev. 108 shows a 
steep rise of variability at E-4.5 keV with a peak at E-5.0 keV 
(Ponti et al. 2004). We analysed the 2 observations with a finest 
time resolution (i.e. 2.5 ks) and found excess variability in all 
the three bands (E-5.4-7.2 keV) during rev. 108B. The energy 
bands we analysed exclude the strong variability event observed 
in the RMS spectrum. 

The variations registered in the entire Fe K band during rev. 
108B are most probably driven by the strong flare observed in 
the continuum light curve (Ponti et al. 2004; Iwasawa & Miniutti 
2004). Papadakis et al. 2005 studied the frequency-resolved 
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Table 4. Variability significances obtained from excess map analysis in the three energy bands: 5.4-6.1 keV (band A), 6.1-6.8 keV 
(band B) and 6.8-7.2 keV (band C). Confidence levels >90% are marked in boldface. 



Sources 


Band A 

(%) 


BandB 

(%) 


BandC 

(%) 


Sources 


Band A 

(%) 


BandB 

(%) 


BandC 

(%) 


IC 4329a (rev. 210) 


35.8 


99.7 


71.2 


MRK 590 (rev. 837) 


71.3 


11.6 


6.8 


IC 4329a (rev. 670) 


98.5 


74.5 


11.2 


















MRK 766 (rev. 82) 


20.3 


32.2 


45.0 


MCG -5-23-16 (rev. 363) 


31.0 


77.3 


66.9 


MRK 766 (rev. 265) 


1.4 


80.2 


49.5 


MCG -5-23-16 (rev. 1099) 


63.2 


79.8 


89.5 


MRK 766 (rev. 999) 


70.8 


38.8 


71.0 










MRK 766 (rev. 1000) 


45.6 


60.7 


3.1 


NGC 3783 (rev. 193) 


36.7 


25.6 


82.7 


MRK 766 (rev. 1001) 


47.6 


59.8 


70.8 


NGC3783 (rev. 371) 


68.1 


5.3 


58.3 


MRK 766 (rev. 1002) 


84.7 


70.0 


66.1 


NGC 3783 (rev. 372) 


96.7 


5.3 


65.2 


MRK 766 (rev. 1003) 


97.2 


84.6 


29.9 










MRK 766 (rev. 1004) 


95.7 


95.6 


24.3 


NGC 5548 (rev. 191) 


61.9 


6.8 


94.3 










NGC 5548 (rev. 290) 


23.0 


82.9 


74.7 


NGC 7469 (rev. 192 A) 


42.9 


4.6 


96.6 


NGC 5548 (rev. 291) 


77.3 


51.1 


14.9 


NGC 7469 (rev. 192B) 


81.6 


40.5 


94.5 










NGC 7469 (rev. 912) 


84.5 


86.9 


86.7 


NGC 35 16 (rev. 245) 


30.2 


28.6 


98.2 


NGC 7469 (rev. 913) 


66.6 


43.3 


3.4 


NGC 3516 (rev. 352) 


88.6 


29.6 


6.5 










NGC 35 16 (rev. 1250) 


14.4 


88.7 


95.7 


MCG -2-58-22 (rev. 180) 


54.4 


44.9 


49.4 


NGC 35 16 (rev. 1251) 


91.7 


62.8 


93.7 










NGC 3516 (rev. 1252) 


67.9 


39.8 


99 


NGC 526A (rev. 647) 


78.7 


58.7 


52.3 


NGC 3516 (rev. 1253) 


77.0 


4.2 


66.6 


















inoi^ wji ^rev. zo j) 


65 4 


57 9 


78 3 


MRK509 (rev. 161) 


95.3 


90.4 


42.3 


NGC 4051 (rev. 541) 


20.3 


53.4 


82.3 




66 1 


3 1 


72 2 










MRK 509 (rev. 1073) 


92.7 


13.7 


88.1 


ESO 141-G055 (rev. 1435 A) 


93.4 


9.0 


93.0 


MRK 509 (rev. 1074) 


83.7 


30.0 


25.1 


ESO 141-G055 (rev. 1435 B) 


13.0 


12.8 


20.8 


iviivrv jyjy ^icv. i loo ) 


27 1 


91.2 


65 8 




19 1 


20 8 


15 8 










ESO 141-G055 (rev. 1445) 


60.7 


59.8 


10.7 


ivii—o -o-ju-ij ^rev. juij 


97.0 


82 5 


93.7 










MCG -6-30-15 (rev. 302) 


63.7 


96.5 


82.8 


UGC 3973 (rev. 1247) 


70.6 


37.8 


52.3 


MCG -6-30-15 (rev. 303) 


83.1 


67.4 


4.7 


UGC 3973 (rev. 1263) 


94.6 


41.5 


74.7 


MCG -6-30-15 (rev. 108A) 


81.6 


79.6 


65.4 


UGC 3973 (rev. 1332) 


37.3 


48.4 


16.2 


U JU lJ ^ICV. 1UOD ) 


97.0 


96.4 


94.5 


uul j?/ j ^rev. ijjjj 


55 6 


24 7 


95.9 


Mrn -1-8 11 11 (rpv 1Q£\ 


14 1 


29 9 


20 3 


Fpirall Q (rt±\T 1 OS"> 


71 8 


25 6 


82 9 


NGC 7314 (rev. 256) 


96.2 


43.3 


67.2 


ESO 198-G024 (rev. 207) 


92.4 


54.7 


35.4 


NdC 71 1 & (rp\r 1 1 7? > 

IN UV^ / Jit ^ICV. 1 1 / A) 


70 9 


115 


17 1 




86 7 


27 6 


2 7 


AKN 120 (rev. 679) 


62.5 


44.0 


45.9 


MRK 110 (rev. 904) 


27.7 


5.6 


56.8 


MRK 279 (rev. 1087) 


64.1 


44.0 


66.0 


AKN 564 (rev. 930) 


12.1 


94.0 


95.5 


MRK 279 (rev. 1088) 


17.5 


69.1 


31.1 










MRK 279 (rev. 1089) 


14.1 


3.6 


94.5 


NGC 7213 (rev. 269) 


64.7 


79.2 


67.1 


MR 2251-178 (rev. 447) 


58.2 


72.3 


83.4 


ESO511-G030 (rev. 1042) 


44.3 


87.5 


11.8 


NGC 3227 (rev. 178) 


72.4 


26.4 


32.2 


MRK 704 (rev. 1074) 


85.7 


91.3 


15.5 


NGC 3227 (rev. 1279) 


40.9 


53.4 


17.4 


MRK 704 (rev. 1630) 


36.8 


58.2 


78.0 


IRAS 05078+1626 (rev. 1410) 


24.3 


21.7 


9.8 


NGC 4593 (rev. 465) 


9.2 


9.9 


54.0 



spectra (Revnivtsev et al. 1999) of MCG -6-30-15 during rev. 
108-303 and found that, although the broad Fe Ka line shows no 
significant variations on time scales less than -1-2 days, some 
residuals between 5-7 keV are suggestive of small amplitude 
variations in a line-like feature. Moreover the drop observed in 
the spectra in the Fe Ka regime is less than the one expected 
from a constant line. Vaughan & Fabian 2004, using the flux-flux 
linear correlation (between soft band, 1-2 keV, and hard band, 
3-10 keV) technique, deduced that flux variations are dominated 
by changes in the normalization of the power law component. 
This is in agreement with our assumption in the Monte Carlo 
simulations, where we assume that only the power law normal- 
ization varies, while the other components are constant. In ad- 
dition there must be an additional component that varies little, 
and contributes more in the 3-10 keV band than in the 1-2 keV. 
The largest fractional contribution from this component occurs 
in the Fe K band, where it accounts for -40% of the total flux. 
Once the variable continuum is subtracted from the total spec- 
trum (see Sect. l3.3D , the residuals should include only the nearly 



constant reflection component. Hence, the significant variations 
we observed in these residuals (during rev. 108B, 301 and 302) 
must be due to additional line-like features overlapping the broad 
Fe line profile, in agreement with deductions of Papadakis et al. 
2005. 

NGC 7314: a relativistic line or a distant reflector model 
can fit equally well the spectrum of this source during rev. 256 
(NOGR07). For our analysis we assumed a simple model which 
accounts for the neutral and highly ionized Fe Ka emission 
lines. We found significant variations redward the neutral com- 
ponent during rev. 256. The detection of a narrow emission fea- 
ture at E-5.8 keV was also claimed by Yaqoob et al. (2003), in 
a Chandra observation taken 19-20 July 2002, and interpreted 
as redshifted Fe Ka emission. A hot spot origin of the line was 
ruled out by the authors because the estimated spot location re- 
sulted too large (i.e. r>132 r^) to cause sufficient gravitational 
and Doppler redshift. In our spectrum, however, the feature is 
observed for about 20 ks (i.e. about half the duration of the 
observation). If produced in a spot - assuming a SMBH mass 
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of 5xl0 6 M o (Padovani & Rafanelli 1988) - the emitting radius 
should be r>26 in order for the emission to correspond to less 
than one orbit (i.e. a spot and not an annulus). This estimated 
lower limit on the radius is still within the upper limit of 132 r g . 
It is worth noting that the target of the second XMM-Newton 
observation analysed in this paper (rev. 1172) is the cluster 
XMMUJ2235.3-2557, while NGC 7314 is observed off-axis, re- 
ducing the X-ray telescope's effective area. 

NGC 3227: We detected a redshifted narrow emission fea- 
ture during the longest observation (rev. 1279). Estimated pa- 
rameters of this feature are in agreement with those derived in 
Markowitz et al. 2009 and double checked with MOS data. The 
energy of the feature results inconsistent with a Compton shoul- 
der minimum energy (i.e. E>6.24 keV). We did not find any sig- 
nificant evidence for variability in the Fe K band, in agreement 
with results of time resolved analysis in Markowitz et al. 2009. 

MKN 766: A relativistic Fe K line was detected by NGO07 
during rev. 82 and 265. To fit the curvature induced by this line 
we adopted a broad Gaussian component whose centroid en- 
ergy falls in the energy range typical of a highly ionized Fe Ka 
and/or Fe KJ3 line. Turner et al. 2006 reported on significant en- 
ergy modulation (detected via excess map technique) during rev. 
265, within E~6. 1-6.8 keV (band B), tracing a sinusoid of pe- 
riod -150 ks. They associated this emission to orbiting material 
at an estimated distance of r=115 r^ from the central BH. At 
such distance the increase in flux during the approaching phase 
of the orbital motion is not expected to be strong (e.g. Dovciak 
et al. 2004), in agreement with the results of our temporal anal- 
ysis, which rules out the presence of any significant intensity 
modulation in the data. An extensive analysis of the entire set of 
observations (rev. 82, 265, 999-1004) available for this source, 
was carried out by Miller et al. 2006. They found variability in 
the low-ionization band (E=6. 08-6.48 keV), significant at 95%. 
This value is in agreement with the variability significance we 
estimated in Band B (which includes the low-ionization band of 
Miller et al.) during rev. 1004. 

NGC 4051: We found evidence for broad emission at Fe 
K redshifted energies in both the observations (rev. 263 and 
541). For what concerns rev. 263, our result is consistent with 
NOGR07. We did not detect significant variability in this source, 
in agreement with the RMS spectrum (time resolution 2 ks) anal- 
ysis carried out by Ponti et al. 2006. 

AKN 564: ASCA observations of AKN 564 revealed an Fe 
Ka line origin in highly ionized material (dominated by H-like 
ions). The line results significantly broad (Turner et al. 2001). 
Variability measurements down to time scales of approximately 
a week constrained the bulk of the line to originate at small dis- 
tances from the nucleus. Previous studies (Papadakis et al. 2007) 
on the -100 ks XMM-Newton observation (rev. 930) revealed a 
Fe Ka line component slightly broader than the energy resolu- 
tion of the EPIC-pn detector and rather weak (i.e. EW-70 eV), 
in agreement with our best fit model (Table [3]). However the line 
cannot be unambiguously ascribed to relativistic emission (see 
Papadakis et al. 2007; NOGR07). We detected significant vari- 
ability in both band B and C where the only residuals visible in 
the map coincide with the Fe line. They seem to be absent at the 
very beginning of the observation; we rule out this to be a back- 
ground subtraction effect, because the source signal results very 
high with respect to the background level (during the first 10 ks 
the background count rate in the 4-9 keV energy band is just 5% 
the source one). 
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